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Abstract. Presented are complete electroweak (EW) corrections at O(aa^) to top-squark pair pro- 
duction at the Large Hadron Collider (LHC) within the framework of the Minimal Supersymmetric 
Standard Model (MSSM). At this order, effects from the interference of EW and QCD contributions 
have to be taken into account. Also photon-induced top-squark production is considered as addi- 
tionnal partonic channel which arises from the non-zero photon density in the proton. Furthermore, 
the impact of MSSM parameters on the EW corrections is analyzed. 

PACS. I2.I5.Lk Electroweak radiative corrections - 13.87.Ce Production 



1 Introduction 

The lighter top-squark is as a candidate for the lightest 
squark within many supersymmetric models [1 , for 
mainly two reasons based on the large top Yukawa 
coupling. Evolving the scalar masses from the GUT 
scale to low scales leads to a low value of the top- 
squark mass, and, moreover, a large mixing in the top- 
squark sector induces a substantial splitting between 
the two mass eigenstates [2^ . Top squarks are therefore 
of particluar interest, especially for hadron colliders. 

In hadronic collisions, top-squarks are primarily 
produced in pairs via the strong interaction. Present 
experimental limits from the Tevatron RUN II data, 
dependent on the lightest neutralino mass, are pro- 
vided by the CDF and D0 collaborations [3]. Con- 
cerning the theoretical predictions, Born- level cross 
sections calculated in [3], have been improved by in- 
cluding the next-to-leading (NLO) corrections in su- 
persymmetric QCD (SUSY-QCD). These were worked 
out in [5] with the restriction to final state squarks of 
the first two generations. The analysis for the stop 
sector, performed in [6], shows that the SUSY-QCD 
corrections significantly modify the LO cross section. 

At lowest order in QCD as well as at 0(a^), only 
diagonal top-squarks pairs can be produced. The non- 
diagonal production is suppressed as the cross sec- 
tion becomes non-zero only at 0(af). The production 
of non-diagonal top-squark pairs can also proceed at 
0(a 2 ) via Z-exchange in e + e~ annihilation [7], or qq 
annihilation [8] . The LO cross section for the diagonal 
pair production depends only on the mass of the pro- 
duced squarks. As a consequence, bounds on the cross 
section can easily be translated into lower bounds on 
the lightest top-squark mass. At NLO, the cross sec- 
tion is not only considerably changed, but also other 
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supersymmetric parameters, like mixing angle, gluino 
mass and other sparticle masses, enter through higher 
order corrections. On the other hand, once the top- 
squarks are discovered, their masses could be directly 
determined from the cross section measurement. 

In the following, we study the NLO EW-like cor- 
rections to the top-squark pair production within the 
Minimal Supersymmetric Standard Model (MSSM). 
We assume the MSSM with real parameters, R-parity 
conservation and minimal flavor violation. 



2 EW NLO Contributions 

As a subset of the complete set of EW virtual correc- 
tions, photonic contributions are present. These con- 
tain IR singularities, which cancel when also the real 
photonic corrections are taken into account. In addi- 
tion, at NLO a photon-induced subclass of contribu- 
tions appears as an independent production channel. 

2.1 Virtual Corrections 

The virtual corrections can be classified according to 
self-energy, vertex, box, and counter-term contribu- 
tions dressing the Born-level partonic amplitudes of 
qq annihilation and gluon fusion. Getting an UV finite 
result requires renormalization of the involved quarks 
and top-squarks. The counterterms for self-energies, 
quark and squark vertices and the squark quartic in- 
teraction at one-loop order are determined in the on- 
shell renormalization scheme. It is not necessary to 
renormalize the gluon field and the strong coupling 
constant. 

Loop diagrams involving virtual photons generate 
IR singularities. According to Bloch-Nordsieck [jJJ, IR- 
singular terms cancel against their counterparts in the 
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Fig. 1. Feynman diagrams for gluon bremsstrahlung con- 
tributions (left) that have to be combined with EW Born 
diagrams (right). Only interference between initial and fi- 
nal state gluon radiation is non- vanishing. 



real photon corrections. To regularize the IR singular- 
ities we introduce a fictitious photon mass A. In case 
of external light quarks, also collinear singularities oc- 
cur if a photon is radiated off a massless quark in the 
collinear limit. We therefore keep non-zero initial-state 
quark masses m q in the loop integrals, which give rise 
to single and double logarithmic contributions of quark 
masses. The double logarithms cancel in the sum of 
virtual and real corrections, single logarithms, how- 
ever, survive and have to be treated by means of the 
factorization. 

An additionnal source of IR singularities originates 
from the gluonic insertions to the box contributions 
in the qq channel. These appear in combination with 
photons or Z-bosons. The similiarity between gluon 
and photon in the box contributions allows us to treat 
the gluonic IR singularities in analogy to the photon 
case. 



2.2 Real Corrections 

To compensate IR singularities in the virtual EW cor- 
rections, real photonic and gluonic contributions are 
required. In case of gg fusion, only photon brems- 
strahlung is needed, whereas in the qq annihilation 
channel, also gluon bremsstrahlung of the appropriate 
order of O(aaf), shown in Fig.Hl has to be taken into 
account. The necessary contributions originate from 
the interference of QCD and EW tree level diagrams 
which vanishes at LO. Not all of the interference terms 
contribute, however. Owing to the color structure, only 
the interference between initial and final state gluon 
radiation is non-zero. 

Including the EW-QCD interference in the real 
corrections does not yet lead to an IR- finite result. Also 
the IR-singular QCD-mcdiated box corrections inter- 
fering with the 0(a) photon and Z-boson tree-level 
diagrams are needed. Besides the gluonic corrections 
there are also IR-finite QCD-mediatcd box corrections, 
which contain gluinos in the loop. Interfered with the 
0(a) tree- level diagrams, these also give contributions 
of the respective order of 0(aa 2 ). 

We encounter also IR-finite bremsstrahlung contri- 
butions, which are, however, suppressed by more than 
a factor of ten and therefore not included in our numer- 
ical studies. For similar reasons we also neglect contri- 
butions coming from the interference of IR-singular 
and IR-finite terms. 



The treatment of IR-singular bremsstrahlung is done 
using the phase space slicing method. The photonic 
(gluonic) phase space is split into soft and collinear 
parts, which contain singularities and into non-collinear 
hard part, which is free of singularities and can be inte- 
grated numerically. In the singular regions, the squared 
matrix elements for the radiative process factorize into 
lowest-order matrix elements and universal factors con- 
tainig the singularities. 

The soft-photon part of the radiative cross section 
in the qq annihilation channel, which is similar to that 
in e+e~ -> ti [TO]. 
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and for the radiative cross section in the gg fusion 
channel, 
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can be expressed using universal factors, o so j t , 
which refer to initial state radiation, final state ra- 
diation, or interference of initial and final state radia- 
tion, respectively, with e q and e* denoting the electric 
charges of the initial-state quark and of the top-squark, 
respectively, da^' 39 denote the corresponding partonic 
lowest order cross sections. 

The soft-gluon part for the qq channel can be writ- 
ten in a way similar to |T]). but with a different ar- 
rangement of the color matrices, 
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with M. denoting the "Born" matrix elements for g, 7 
and Z exchange where the color matrices are factorized 
off. 

The collinear part of the 2-^3 cross section is 
proportional to the Born cross section of the hard pro- 
cess with reduced momentum of one of the partons. 
Assuming that parton a with momentum p a radiates 
off a photon with p 7 = (1 — z)p a , the parton momen- 
tum available for the hard process is reduced to zp a . 
Accordingly, the partonic energy of the total process 
inclusive photon radiation s — (p a + pb) 2 , is reduced 
for the hard process to s = (zp a + Pb) 2 ■ 

Using these variables, the partonic cross section in 
the collinear cones can be written as [TTl[r2] 



da C ou (s) 



With Kcoll(z) = 7^ P qq( Z ) 



a- 
— 

7T 
1 

2 J 



1-<5 S 



+ 2 (1 - Z) < 



dz da (s) K C oii{z) , 



(4) 













- 1 









where P q q(z) = (1 + z 2 )/(l — z) is an Altarelli-Parisi 
splitting function [13] and 5g is the cut-off parameter to 
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define the collinear region by cos 8 > 1 — 8g ■ The Born 
cross section refers to the hard scale s, whereas in the 
collinear factor the total energy s is the scale needed. 
In order to avoid an overlap with the soft region, the 
upper limit on the z-integration in ^ is reduced from 
z = 1 to z = 1 — S s , where S s — 2AE/y/l. 

As already mentioned, after adding virtual and real 
corrections, the mass singularity in ((4]) does not cancel 
and has to be absorbed into the (anti-)quark density 
functions (PDFs). This can be formally achieved by a 
redefinition at NLO QED as shown in [TTlll4ill5j . 



2.3 Photon-Induced Top-Squark Pair Production 

We also consider the photon-induced mechanism of 
top-squark pair production, which becomes non-zero 
at NLO in QED as a direct consequence of the non-zero 
photon density in the proton. Although the photon- 
induced processes are of different overall order, they 
contribute to the same hadronic final state and thus 
represent contributions at NLO QED. We consider 
only the photon-gluon process in our numerical analy- 
sis and neglect the quark-photon process, which gives 
contribution of higher order. 

As the PDFs at NLO QED have become available 
only recently [16j . the results shown here thus corre- 
spond to the first study of these effects on the top- 
squark pair production [T7] . 



3 Numerical results 

Here we present numerical results for the production of 
lighter top-squark pairs at LHC energies with EW con- 
tributions at one-loop level. We show in the integrated 
hadronic cross section a and the differential hadronic 
cross sections with respect to the invariant mass of the 
top-squark pair inclusive the photon (do~/dMi nv ) and 
to the transverse momentum (da/dpr) of one of the 
final state top-squarks. All hadronic quantities are ob- 
tained by folding the partonic cross section with par- 
ton distributions and summing over all contributing 
partons. 

Our Standard Model input parameters are chosen 
in correspondance with [19]. As described above, we 
use the MRST 2004 QED [TB] PDF set with the choice 
of factorization and renormalization scales to equal the 
sum of the final state particles, /Up = fiR = 2mt. qq 
denotes the sum of the uu, dd, cc, and ss annihilation 
channels. 

In Table Q] we show results for the total hadronic 
cross sections within four different SPS mSUGRA sce- 
narios p~8l[T9] . The total cross sections at leading order, 
a LO , the absolute size of the EW corrections, which 
corresponds to the difference between the LO and NLO 
cross sections, Aa , and the relative corrections, S, 
given as the ratio of NLO corrections to the respec- 
tive LO contributions, are presented for gg fusion, qq 
annihilation, and 177 fusion separately. The 57 chan- 
nel contributes only at NLO. The total amount of EW 
corrections is small in all SUSY scenarios. 



Table 1. Numerical results for the top-squark pair produc- 
tion at the LHC within different SPS mSUGRA scenarios 
18, 19] (the values of a and Aa are given in fb). 



scenario 


channel 


a LO 


Aa NLO 


Aa NLO 


SPSla 


qq 
gg 
y / 
total 


220 
1444 

1664 


-9.65 
-15.4 
29.0 
3.95 


-4.4 % 
-1.1 % 

0.24% 


SPSla' 


qq 
gg 
y / 
total 


436 
3292 

3728 


-11.5 
-14.6 
58.5 
32.4 


-2.6 % 
-0.44% 

0.87% 


SPS2 


qq 
gg 

91 
total 


1.16 
2.97 

4.13 


-8.98xl0 -2 
-3.07xl0" 2 
0.155 

3.45xl0 -2 


-7.7 % 
-1.0 % 

0.84% 


SPS5 


qq 
gg 

91 
total 


2870 
31960 

34830 


-13.2 
499 
405 
891 


0.46% 
1.6 % 

2.6 % 
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Fig. 2. Overall relative corrections Stot with respect to pr 
and ■s/S at the LHC within the SPSla scenario. 



In order to illustrate the numerical impact of the 
EW corrections on the LO cross section, in Fig. [2] we 
show the overall relative corrections Stot = Aa NLO / a LO 
as distributions with respect to px and y/l (which cor- 
responds to M inv ). In the pj-distribution, the correc- 
tions grow in size with increasing pt and reach about 
-20% for p T > 2500 GeV. Similar effects are visible 
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in the Mi nv distribution. Although smaller in size, the 
corrections raise up to —15% level for \f§ > 5000 GeV. 
These ranges of pr and y/l are still within the reach 
of the LHC. 

The behavior of EW corrections at high scales is 
dominated by the massive gauge boson contributions, 
which consist of double logarithms of W and Z masses. 
The double logarithms are not canceled by the real 
gauge boson corrections, since these correspond to dif- 
ferent hadronic final states. As a result, large negative 
contributions show up in the pt and y/§ distributions. 

It is obvious that although small for the total cross 
sections, the EW corrections cannot be neglected in 
the differential distributions, since in the high-p^ and 
high-\/i ranges they are of the same order of magni- 
tude as the SUSY-QCD corrections. 

We have also studied the dependence of the EW 
contributions on various SUSY parameters. We have 
focused on the parameters that affect the top-squark 
mass, since here the effects are expected to be strongest. 
Following parameters have been varied: tIq 3 , m (j 3 > 
tan/3, A t , and fi around the SPS la' value, while keep- 
ing all other parameters fixed. As an example, in Fig. [3] 
we show the dependence of the overall EW corrections 
on ?tiq 3 for each production channel separately. 

The 37 corrections grow up to 2% with increasing 
mj in the considered range. The gj fusion channel is 
as important as the qq and gg channels and should be 
taken into account for a reliable cross section predic- 
tion. The qq corrections involve many different SUSY 
particles in the loops, however the relative corrections 
change only little between different scenarios, varying 
between 0% and —1%. The gg contributions are more 
sensible to variations of the considered SUSY param- 
eters. The general behavior is similar to the qq case, 
the decrease with increasing mi is, however, stronger. 
The gg plot is dominated by striking (negative) peaks, 
some of them are (although much weaker) also visible 
in the qq corrections. 

The peaks originate from two sources of threshold 
effects. One of them is the Higgs boson H° threshold, 
m(H ) — 2mg , which affects only the gg channel. Sec- 
ond source of the threshold effects enhances the correc- 
tions in scenarios where equals the sum of masses 
of a neutralino and the top-quark or of a chargino 
and the bottom-quark. The parameter regions with 
large EW corrections due to the threshold effects are 
small and over the wide range of SUSY parameters, 
the EW corrections to the top-squark pair production 
are smaller than 1%. 
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